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^ We evaluate the interaction potential between a hydrogen and an antihydrogen using the second- 

order perturbation theory within the framework of the four-body system in a separable two-body 
basis. We find that the H-H interaction potential possesses the peculiar features of a shallow local 
minimum located around interatomic separations of r ~ 6 a.u. and a barrier rising at r < 5 a.u. 
Additional theoretical and experimental investigations on the nature of these peculiar features will 
be of great interest. 

Q ' PACS numbers: 31.15.Md, 31.15.-p, 34.10.-l-x, 36.10.-k, 

• i-H ' 

I. INTRODUCTION 

^ ' Recent production of cold antihydrogen [1, 2] has stimulated new theoretical and experimental interests on the 
\ H and the H system. The four particles of the greater (e+e~pp) system can be arranged in different ways leading 
. to different types of states. On the one hand, there are states in the H-H family in which the e"*" and the e~ are 
\l ■ not correlated but are orbiting around the p and the p respectively, with the H and the H appearing as composite 
particles. They are the lowest-energy states at large interatomic separations. On the other hand, there are states of 
OO the Pspp family for which the e"*" and the e~ are highly correlated and orbit each other as in a positronium (Ps). 

The role of the Pspp family and the H-H family is interchanged at short distances as the e+-e~-correlated state of 
the Pspp family becomes the state of the lowest energy, as indicated by the merging of the V^ar {f) potential obtained 
by variational calculations with the Vpspp(r) potential obtained for the Pspp configuration, at small interatomic 
separations (see Fig. 4 below). Because the p and p can orbit around (and annihilate) each other under their mutual 
Coulomb and nuclear interactions [3] , the Pspp family can be further bifurcated into subfamilies with a correlated or 
uncorrelated p-p pair. There are also states that are hybrid of different families. 
' If one can depict the orbiting of one particle relative to another particle as a "dance pattern" , then the dance patterns 
of the different families have distinctly different topological structures and connectivities. The distinct topological 
structures may allow them to retain some of their characteristics and stability. The mixing and the interplay between 
different families at different interatomic separations will provide a general idea on the "genealogy" of the states. To 
be able to classify the greater (e+e~pp) system into families, if at all possible, will bring us to a better understanding 
of the complexity of the spectrum that is associated with the complicated four-body problem, as well as to a better 
knowledge of the scattering between the H and the H. 

The H atom and the H atom can scatter from each other at various energies. The problem of the H-H scattering 
will require the knowledge of the interaction potential between the H and the H as composite particles. As the H 
and H may rearrange into other configurations and may annihilate at short distances, the scattering problem needs 
to be treated in general in a coupled-channel analysis involving different interaction potentials and transition matrix 
elements [4]. 

For the elastic channel, one often takes the interaction potential to be the adiabatic potential obtained in a vari- 
ational calculation for the lowest-energy state of the e"*" and the in the Born-Oppcnheimer potential [4-10]. The 
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e'^-er correlation property of the lowest-energy adiabatic state will undergo a change as the interatomic separation 
changers. At large interatomic separations, the adiabatic lowest-energy state coincides with the incoming channel 
state of separated H and H atoms with uncorrelated e~ and e"*" orbiting around p and p, respectively. At small 
interatomic separations, the variational calculations will naturally lead to the lowest-energy adiabatic state involving 
strong correlations between the e+ and the which differs from the elastic incoming channel by the presence of 
e+-e~ correlations. An adiabatic approximation involves altering the nature of the underlying state from the incom- 
ing elastic channel with no e+-e~ correlations to the other e+-e~-correlated channel, as the interatomic separation 
decreases. A general consideration of the channel coupling however will contain both the incoming elastic channel and 
the e"'"-e~-correlated channel in explicit coupling and transition. Depending on the collision energy, the dynamics of 
the scattering process will lead to the adiabatic case with strong e^-e" correlations at short separations in one limit 
of slow motion and large transition matrix elements. It will lead to the diabatic case with no e^-e~ correlation of the 
atomic H and H in the other limit. To carry out the explicit channel coupling for the general case of an arbitrary 
energy, or to see the cross-over to either the; adiabatic or the diabatic limit, it is necessary to obtain the interaction 
potential for the incoming entrance channel oi a, H and an H. 

In the present work, we focus our attention on the interaction potential for this incoming entrance channel between 
the H and the H as composite particles, for which the e~ and the remain uncorrelated, orbiting aroimd the p 
and p respectively, i.e., without the loss of the H and the H atomic identities. We are in effect trying to construct 
the interaction potential for one of the diabatic basis configurations in preparation for a general coupling of different 
diabatic basis states in a couple-channel analysis. For brevity of nomenclature, we shall use the term "the interaction 
potential" to refer to this potential between the composite particles H and H, unless otherwise stated. 

The interaction potential Vj^f^ (r) between the composite H and the H can also be used to investigate possible H-H 
two-body molecular states of the H-H family, as the composite H and H atoms can form a bound system if their 
interaction potential is sufficiently attractive. These possible molecular two-body H-H states, if they exists, belongs 
to the H-H family, a subset of the states of the greater four-body {e^e~pp) system. One theoretical approach to 
study the lowest-energy state of the greater {e'^e~pp) system is to use the variational method equipped with either 
explicitly correlated Slater-type or Gaussian-type orbitals with e+-e~ correlations. This method has provided a great 
wealth of information regarding the lowest-energy state of the four-body system [5 9]. As our interest is focused on 
the properties of possible two-body states of the composite H and H atoms, the two-body H-H states of our interest 
may or may not necessarily be the lowest-energy state of the greater (e+e~pp) system. Its relative position depends 
on the interatomic separation. We mentioned earlier that at large interatomic separations, the two-body state with 
separated atoms and leptons in the Is orbitals is the lowest-energy state of the greater (e+e~pp) four-body system. 
At small interatomic separations, the two-body molecular state of the H and H lies higher than the variational 
adiabatic lowest-energy state which is characterized by strong e^-e~ correlations, and the two-body molecular state 
of the H and H is an excited state of the greater {e^e^pp) system. It will be of great interest in future work to 
investigate how the two-body H-H state may preserve or otherwise mix their distinct molecular characteristics with 
the e"'"-e~-correlated state, as the system dynamically traverses from large interatomic separations to the region of 
small separations. 

To follow the four-body problem in the two-body atomic basis, we study the system from outside-in as the H 
and the H approach each other, and we examine their virtual excitations into excited H-H configurations, as the 
interatomic distance decreases. A good method to carry out such an investigation is to use the perturbation theory, 
having the atomic H and H as unperturbed composite particles, as in previous study of molecular states in heavy 
mesons in hadron physics [11]. The perturbation theory can be justified as a useful tool here because the composite 
H and H atoms arc neutral objects, and their residue interaction Vj between the H and the H contains four terms 
(as given in Eq. (3) below) involving many cancellations due to the opposite charges of e~ and p of onc^ atom on the 
one hand, and additional cancellations due to the opposite charges of the e+ and p of the other antiatom on the other 
hand. As a consequence, the residue interaction is small in magnitude. Using Vi/{AE), where Vj = rg. ■ Vb/r^ is the 
dipole-dipole interaction and AE is the energy denominator from the ground to the first excited state of the H atom, 
we roughly estimated Vi/AE at interatomic separation r = 3, 5 and 8 a.u. to be ~ 0.05, 0.01 and 0.003, respectively. 
Clearly, the next higher-order term involving the ratio (Vi/AE)'^ will be even smaller and the perturbation series 
is expected to converge. Therefore these numbers justify the application of the perturbation theory to the H-H 
interaction potential in intermediate interatomic distances (i.e., 3 < r < 10 a.u.). 

Within the perturbation theory, the H-H interaction potential is at present known only within the leading order 
[12]. It is positive for large distances, with a barrier whose peak lies in the region of short distances. It decreases 
precipitously at very short distances. Previous analysis of the next-to- leading order interaction potential between the 
H and the H assumed an expansion of the energy denominator up to the first order of the state energies [12]. Such 
an expansion may not be accurate. It is of interest to re-evaluate more accurately this interaction potential, Vjjf[{r), 
up to the next-to-leading order (second order) in the residue interaction Vj. We shall be interested in the interaction 
potential for the elastic channel for which the initial incoming and the final outgoing H and H atoms are in their 
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respective ground states at asymptotic separations. Interaction potential for states with excited H and H atoms at 
asymptotic separations can be similarly considered in a simple generalization. 

For small interatomic distances, the e^-e~ correlations become considerably important and we shall also study this 
additional e^-e~ correlation effects on the potential energy around this region of separation. Needless to say, it is 
necessary in the future to include a greater number of the degrees of freedom, to take into account addition types of 
polarization of the two-body states. The perturbation theory results obtained here can provide a benchmark against 
which additional distortions may be measured. They can bring us to a better understanding of the scattering potential 
and the complexity of the spectrum that is associated with the four-body system. 

This paper is organized as follows. In Section II, we review the formulation of the four-body problem in terms of 
the interaction of two composite objects using separable two-body basis. In Section III, we show how to evaluate 
the interaction matrix elements and the interaction potential. The resiilts of the H-H interaction potential and the 
examination of the stability with respect to e+-e^ correlations are given in Sections IV and V, respectively. Finally, 
Section VI gives some discussions and conclusions of the present work. Some of the details of the analytical formulas 
are presented in the Appendix. Atomic units m=h=e=l are used throughout the paper unless otherwise indicated. 



II. THE FOUR-BODY PROBLEM IN A SEPARABLE TWO-BODY BASIS 



We shall review the formulation of the four-body problem in a separable two-body basis in terms of the interaction 
of two composite particles as presented previously [11]. As applied to the H-H system, we choose the four-body 
coordinate system as shown in Fig. 1 and label constituents p, e~, e+, and p as particles 1, 2, 3, and 4, respectively 
with a non-relativistic Hamiltonian 

j=i J j=i k>j j=i 

in which particle j has a momentum pj and a rest mass mj. The pairwise interaction Vjk{rjk) between particle j and 
particle k depends on the relative coordinate between them, rjk = Vj — Vk- We introduce the two-body momentum 




FIG. 1: The coordinates of the {p, e , and p} system. 



Pjk = Pj+Pkj ^-iid the two-body internal relative momentum p^-^, = fkPj — fjPk^ where fk = mk/mjk, rrijk = rrij+mk- 
We partition the Hamiltonian so that in the lowest order the state in question can be described by a state of the 
unperturbed Hamiltonian. We choose to partition H into an unperturbed Hamiltonian hi2 + hzi of atoms ^(12) and 
5(34) according to 

H = + + Vi + hi2 + hs4, (2) 

2mi2 2m34 



Vi = yi4(ri4) + Vi3{ri3) + V23{r23) + V24{r24), (3) 

'ijk = ^+Vjk{rj-rk) + mjk for (jfc) = A(12) and 5(34), (4) 
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where jjjk = mjmk/nijk- The eigenvalues of the Hamihonians hi2 and h^A can be solved separately to obtain the 
bound state wave functions and masses Mjk{y) of atoms A and B, 

hjk\{jk)u) = [ejkii^) +'mjk]\{jk)^) = Mjk{i/)\{jk)„). (5) 

The four-body Hamiltonian becomes 

H=P^ + P^ + Vi + M,2{u)+M34iy)- (6) 
2rni2 2rn34 

The above non-relativistic approximation is obtained from relativistic results by neglecting terms of order tj^/Mj^ 
and higher [13]. In order to satisfy the boundary condition at large separations for which Vi approaches zero, we need 
to include somc^ of these higher-order terms and modiiy mi2 and 77134 in the above equation to Mi2(z^) and M3i{y') 
so that the Hamiltonian becomes 

The above Hamiltonian then describes properly the asymptotic behavior at large separations of the two atoms. This 
Hamiltonian H is related to the standard Hamiltonian 

= £ + ^ + E E + E (8) 



2me 2me+ 

hy H = iJstand + <^12{v) + €34(1^') + A, 



0=1 k>j j=l 



A = ^ + ^ + ( — \ i ^ — \ (9) 

2mp 2mp \2M12i1y) 2m, 12 J 2^34(1^') 2m34 / ' ^' 

and A is a small quantity which can be neglected in the Born-Oppenheimer limit of fixed baryon centers. 

We consider a four-body system in which both atoms will be in their ground states |^o(12)-Bo(34)) at asymptotically 
large interatomic separations where the residue interactions Vj vanishes. When we include Vj as a perturbation, the 
eigenfunction of H becomes 

,T,. ^ u DN , \AxBy){AxBx'\Vi\AoBo) \ 

" ^^^^ = 1 1^°^°^ + g e4A)+e.(A0-e.(O)-e.(O) | ' ^''^ 

where r = R12 — R34 is the interatomic separation (sec Fig. 1). Rjk = fj Vj + fk Vk is the center-of-mass coordinate 
of rrij and rufe, and YL'w indicates that the sum is over all atomic states, including continuum states, except \AqBq). 
The eigenvalue equation is 

F4'(r,ri2,r34) = [Mi2(0) + M34(0) + e]«'(r, ri2, r34). (11) 

Working in the center-of-mass frame and taking the scalar product of the above equation with \AqBq), we obtain the 
Schrodinger equation for the motion of ylo(12) relative to i3o(34), 

^ +y(r))v(r-)=eV;(r), (12) 
where p is the relative momentum of the composite particles 

^ M34(0)Pl2-Mi2(0)P34 
^ Mi2(0) + M34(0) ' ^ ' 

and IJ.AB is the reduced mass of the two atoms 

_ Mi2(0)M34(0) 

- Mi2(0) + M34(0)- ('^^ 
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III. EVALUATION OF H AND H INTERACTION MATRIX ELEMENTS 



The Schrodingcr equation (12) involves the interatomic separation r and describes the motion between two com- 
posite atoms. Because the mass of the proton is much greater than the mass of the electron and the positron, the 
interatomic separation r is approximately equal to the internuclear separation ri4. The atom-atom potential V{r) 
in Eq. (12) is given by 

^A,A' e^(A) + eB(A') - e^(0) - 6^(0) 

where \AqBq) is a product of two Is orbitals of the hydrogen and the antihydrogen. We label V{r) as Vjjfj{r), and 
introduce E{r) = V{r) — 1. We call the first (leading-order) term on the right hand side the direct potential, Vdir{r), 

Vdirir) = {AoBo\Vi\AoBo), (16) 

and the second (next-to-leading order) term the polarization potential, Vpoi{r), 

V \{AxBy\Vi\AoBo)\^ 

The direct potential is given by the sum of the four matrix elements {AoBo\Vjk{rjk)\AoBo) of Eq. (3). Prom the 
Is wave functions we can analytically determine [12] 

{AoBo\Vi\AoBo) - ^ (-1 - + + ^^') ■ (18) 

The polarization potential Vpoi (r) is obtained as a double summation of A and A', representing the virtual excitation 
of the H and H from Is states to excited A- A' atomic states. It contains contributions where A- A' are bound-bound, 
bound-continuum (or vice- versa), or continuum-continuum states. 

For the evaluation of the bound-bound component of the polarization potential, Vpl^i{r), we shall use the Fourier 
transform method [11, 14]. We represent the internal wave functions in ^^(12) and i?A'(34) by normalized eigenfunc- 
tions of the hydrogen atom 0f (ri2) and (f)y{r34), respectively. The matrix elements {A\B\/\Vi\AoBo) is then a sum 
of four terms with Vj given by Eq. (3). Each of these terms is given by 

{AxBx,\VMrjk)\AoBo) = J ^e'P-^p^o[/A(ifc)p]pf,o[/B(ifc)p]t;,-,(p), (19) 

where 

pi^{p)= Jdye'Pypi^iy), (20) 

Vjk{p) = J drjke-'P-''='Vjkirjk), (21) 

with p"^Q{ri2)=(j)'{{ri2)4>o{i'i2), Px'o{i'34,)=4'*y{r34)(j)o{r34), and f{A.B}ijk) are the coefficients of the linear relation 

rjk = r + fA{jk) ri2 + fsUk) ru, (22) 

where the values of f{A.B}{.jk) arc given in Ref. [11]. 

To calculate the matrix element (19), we fix the vector r to lie along the z-axis, and quantize the azimuthal 
component of the magnetic quantum number m to be projections along the z-axis. We limit our consideration to 
intermediate states with 1=1 orbital angular momentum, as the dominant polarizing interaction is the dipole-dipole 
excitation. The bound-bound intermediate states are then (A, A') = ({n, (/=1), to}, {n', {1=1), m'}) where n and n' are 
principal quantum numbers and n, n' > 2. The matrix clement {A\Bx'\Vjk{rjk)\A()Bo) is non-vanishing only for 
{A = Im; A' = 1 (—to)} with {to = —1, 0, 1}, when both atoms are excited. 

We write down the Is and np hydrogen wave functions explicitly and obtain the Fourier transform of the transition 
density p^l^ {p) (see Eq.(A9) in Appendix A) 



and 



P.0 iP) - 2^^[ k ) - ^) ^-'^ [^l) ' (23) 
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where /3 = l/(n + 1), 



"2(n+l)" 



(n-2-fc)! 
2n{n+l)\ ' 



(24) 



( = V (n+1)! Ik^T+IiO'^"-^-^ I (#FTlj ^''^ 



Upon substituting the above Fourier transform of the transition density Pxo^{p) into Eq. (19), we can carry out 



and Cl{q) denotes the Gcgcnbaucr polynomials 
Upon substituting the above Fourier transfo: 
the angular part of the integral analytically. We recognize 

oo 

^ipr ^ ^(2/ + iyji{pr)Pi{cosep), (26) 

(=0 

and we can easily carry out the angular integral by using 

27r 1 . 

[ riA [ A„ jp r V /a \v fa \ J ^o{pr) - '^n{pr) (for m = 0), 
J J \3o[pn+32{w) (form = l), 

0-1 

where ji{x) and Pi{z) are a spherical Bessel function and Lcgondre polynomial of a degree of /, respectively. The 
remaining one-dimensional integral in the momentum space can be carried out numerically. 

For the evaluation of the bound-continuum V^^i (r) and continuum-continuum V^^i (r) components of the polarization 
potential, we use the bound and continuum Coulomb wave functions in the configuration space [15] and perform a 
six-dimensional numerical integration involving the residue interaction Vi to obtain {A\B\i\Vi\Af)BQ) . The integrand 
is confined to a finite spatial region as the transition densities Pxo{ri2) and pfo(''34) contain the product of the 
continuum wave function and the ground state wave function. The sum over the continuum states can be carried out 
by discretizing the continuum spectrum which leads to converging results without much difficulty, as the magnitude 
of the matrix element decreases for continuum states with large momenta. We check our results against previous 
calculations of the Cq dispersion-energy coefficient at large interatomic separations. Good agreement of the Cq 
coefficient then allows us to proceed to explore the polarization potential at shorter interatomic separations. 



IV. THE H AND H INTERACTION POTENTIAL 

In the calculation to obtain the bound-bound contribution V^^i (r) to the polarization potential, we include all np 

states up to n < rimax- Fig- 2 shows the bound-bound results of r^Vpgi{r) as a function of the interatomic separation r 
for different Umax number of bound states. Clearly, an increase in rimax for rimax > 10 brings essentially little change 
in the shape but an overall small increase in the magnitude of the polarization potential, indicating that the result is 
sufficiently converged by considering up to n < nmax = 20 states. 

Results for the bound-bound, bound-continuum and continuum-continuum polarization potentials defined in Eq. 
(17) are presented in Fig. 3 as a function of the interatomic separation r. To exhibit the asymptotic behavior of 
Vpol ~ Cq/t^ at large r, the potentials are multiplied by in this figure. The numerical asymptotic values for each 
component are given in the parentheses, to be compared with the values obtained by Eisenschitz and London [16]. 
While r^Vp^^i{r) decreases monotonically as a function of r, the bound-continuum r^Vpl;i{r) and contimimn-continmim 
r^V^^i{r) components have minima at r ~ 6 and r ~ 4 a.u., respectively. The sum of all three contributions leads to 
T^Vpoi{r) decreasing monotonically as a function of r. The result in Fig. 3 also illustrates that the contribution of the 
continuum states to the polarization potential is substantial and cannot be neglected. The degree of the modification 
when continuum intermediate states are included at large r are also known from the perturbation calculation for the 
Cq coefficient of H2 by Eisenschitz and London [16]. Note that the best C% coefficient has been known at least up 
to 8 digits of accuracy. Compare to the old calculation [16], the present Cg = 6.498 yields a better agreement to the 
modern Ce = 6.4990267 dispersion-energy coefficient (i.e., less than ~ 0.05%). 

The functional form of the polarization potential Vpoi (r) in Fig. 3 implies that while the polarization potential varies 
as Ce / at large separations, they vary much more slowly as a function of r at small interatomic separations. This 
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FIG. 2; (a) Plot of Vp^i{r) x as a function of the interatomic separation r for different Umax oil — 1 states, (b) Ratio of the 
summation from 2p up to rimaxP states to the summation from 2p up to 2Qp states as a function of the interatomic separation 




FIG. 3: Bound-bound, bound-continuum and continuum-continuum components of the polarization potential multiplied by r' 
as a function of the interatomic separation r. 



slow variation and the relatively flat behavior can be seen from the total polarization potential Vpoi{r) as a function 
of r shown as the dashed curve in Fig. 4. 

In the present calculation, we preserve the H and H atomic identities and let them interact through Vi defined 
in Eq.(3). Within such a framework, the interatomic potential, Vdir^ Vpoi and the total V^jfj are shown as various 
curves in Figs. 4 and 5. In Fig. 4, the direct potential Vdir from Eq. (18) is slightly repulsive at large distances 
with a barrier whose peak lies in the region of small interatomic separations. It eventually decreases precipitously at 
smaller interatomic separations. The polarization potential Vpoi{r) is, however, always attractive (Fig. 3). The overall 
magnitude of the polarization potential is small because the residue interaction Vj consists of four terms which tend to 
cancel among themselves. The polarization potential provides a small overall modification of the interaction potential 
at short distances (Fig. 4), but it changes the character of the interaction potential at intermediate separations of 
r ~ 6 a.u. where the direct potential is small in magnitude, and the total interaction potential Vfjf[{r) exhibits a 
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shallow local minimum at r ~ 6 a.u. (Fig. 5). 

To understand the rearrangement of the four particles in H and 77, we also plotted in Fig. 4 the potential V^arif) 
based on the variational calculation for the lowest-energy state of the greater (e^e^pp) system [8]. Asymptotically, 
the variational potential Vyarir oo) goes to zero in energy and almost coincides with Vfff[{r) at large interatomic 
separations. At intermediate interatomic separations, the variational potential Vyarir) lies below Vij^f[{r). At very 
short distances (i.e., r < 2 a.u.), Vvar{r) decreases precipitously. 

In Fig. 4, we also present the potential 

Vpsppir) - -— - 0.25 + 1 (28) 

for the Pspp configuration of a positronium and two heavy nuclei. Comparing T4ar(^) with Vpspp{r) at large r, the 
energy gap between them is 0.75 a.u. Only when r is small and the two curves start to meet at r ~ 1 a.u. will 
the rearrangement from H-H to Pspp becomes more probable. In the vicinity of the "critical crossing-distance rc" 
between the Vuar{r) and the Vpspp{r) curves, the Born-Oppenheimer or adiabatic approximation starts to break down 
[9, 17]. Nevertheless, in order to rearrange H-H to become Pspp as r decreases certainly required some changes in the 
topology of the four-body system. At r slightly greater than 1 a.u., the wave function may be regarded as a mixture 
of the two families of HH and Pspp. As the separation decreases further below r < 1 a.u., the Vpspp{r) and Vvar{r) 
potential curves merge, indicating that the Pspp with strong e~^-e~ correlations is the state of the lowest energy. The 
mixing of the two families at small interatomic separations and the strong e+-e^ correlations for r < 1 a.u. is in 
agreement with our analysis of stability against the variation of the e+-e^ correlations in the next Section. 




r (a.u.) 

FIG. 4: The H-H interaction potential V[jf[{r) and its components Vdir{r) and Vpoi{r), as a function of the interatomic 
separation r. Shown here are also the potential Vpspp(r) for the {e'^e~)pp configuration and the potential Vvarir) for the 
lowest-energy {e^e~pp) four-body state obtained from the variational calculations of Ref. [8]. 

In Fig. 5, we show a magnified view of the interaction potential of H-H for interatomic separations from 4 to 10 
a.u. In this region, the sum of the direct potential and the polarization potential gives rise to the H-H interaction 
potential that possess a shallow local minimum located around r ~ 6 a.u. and a barrier rising at r < 5 a.u. The 
depth of the well is found to be 6 x 10~^ a.u. or 1.63 meV, which may hold a bound state. 

It is easy to understand the formation of a potential "pocket" for H-H. The system has the long-range attractive 
potential behaving as ^ Cq/t^ at large r due to the polarization of the atoms as they approach each other. It 
experiences a repulsion at shorter distances arising from the repulsion of the e~ with the p and the 6+ with the p, 
which becomes effective at 1.5 < r < 5 a.u. The combination of a long-range attraction and a shorter range repulsion 
results in a potential pocket at r ~ 6 a.u. The interaction however turns to become attractive at smaller distances of 
r < 1.5 a.u. due to the attraction of the leptons and the baryons. 



9 




r (a.u.) 



FIG. 5: A magnified view of Fig. 4 in the region of intermediate separations. Shown here are the H-H interaction potential of 
Vfjiiir) and its components Vdir(r) and Vpoiif), together with the potential Vvarir) [8]. 

Upon comparing with the variation potential T^arC?"), we observe that Vfjfj{r) almost coincides with Vyar{f) at 
the asymptotic region of large r. As r decreases, Vfj^j deviates from and lies above Vyar{r) of the lowest-energy 
state, indicating that the state of the H-H family is an excited state. At this point, we would like to point out that 
our aim here is to question, when there are no leptonic correlations in our trial wave functions, how the neutral H 
and H interact at intermediate distances before they dissolve completely into positronium and protonium at shorter 
distances. Fig. 5 provides the answer to our question. 

V. STABILITY AGAINST VARIATIONS IN e+ e" CORRELATIONS 

Unlike the interaction between two hydrogen atoms, there is no exchange symmetry between the two leptons in the 
H-H system. However, the correlation between the electron and the positron remains important in the electrostatic 
energy calculations. Therefore it is intriguing to see, within the present approach, when the e^-e~ correlation can play 
an important role in the energy of the four-particle system. It is also crucial to find out whether the local minimum 
at r ^ 6 a.u. examined in the last section is stable against e^-e~ correlation perturbations. We therefore construct 
first a "poor man" variational wave function with the variational parameter ry that allows for e^-e~ correlations, 

^'(r,ri2,r-34,r23) = A^,{y(r^^Ao(ri2)Bo(r34) + ?7Co(r23)i?o(ri4)}, (29) 
where iV^ denotes the normalization constant such that for a fixed r 

j dri2dr34|^'(r,ri2,r34,r23)P = 1 (30) 

and 

j dri2dr34|Co(r23)i5o(ri4)|' = 1. (31) 

The (^0(^23) is the ground-state positronium wave function centered at the center of mass origin and Do('"i4) is a 
plane wave state in the continuum. We shall examine the case of Z?o(''i4) to be a state with zero momentum (for p 
and p to be at rest) which supplies a constant to normalize the product wave function in Eq. (31). Wc numerically 
evaluate (4'|^fstand(f')|^) = E{r) as a six-dimensional integral in the configuration space. 

Figure 6 gives the quantity E[r) = Vuij{r)—1 as a function of the variation parameter rj for different interatomic 
separations r. For the largest values of r in Fig. 6(a), the minimum of E{r) resides at 77 = 0, indicating that the HH 
system is stable against e^-e~ correlations at large r. As r decreases, the E{r) potential surface flattens as a function 
of 77 and the H-H system remains stable against e^-e~ correlations. In the region where the energy minimum begins 
to shift away from r; = 0, the variation of the energy as a function of rj is shown for a finer increment of r in Fig. 6(6). 
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FIG. 6: E{r) as a function of r] and the interatomic separation r. 



At 2.25 > r > 1.2 a.u., the minimum of E{r) is shifted to a positive value of rj. The equihbrium configuration is one in 
which the state is a mixture of the HH and the Pspp configuration. The amphtude r] at equihbrium increases as the 
interatomic separation decreases. As a result, the potential obtained from the variational method, Vvar{i"), is lower 
than Vffjj and Vpspp in this range of separation. At r < 1.2 a.u., the lowest-energy state occurs at = ±1 (Fig. 6(a)) 
which corresponds to a state with a completely correlated e^-e" component Co(r23)£>o(ri4) and the absence of the 
uncorrected component ^o(ri2)-Bo(T"34). The equilibrium configuration resides in the Pspp configuration, indicating 
the dominance of e^-e" correlation at small r. This picture is consistent with the findings of earlier variational 
calculations [5-9] and Fig. 4. 

While Fig. 6 gives the general feature of the potential landscape with respect to e+-e~ correlations as r changes, 
it is necessary to investigate whether the local potential minimum at r = 6 a.u. is stable against perturbations in 
e+-e~ correlations. We replace I^oSq) in Eq. (29) by the wave function represented by the curly bracket of Eq. (10), 
with the \AxBy) amplitudes evaluated at r = 6 a.u. and the summations over A and A' carried out to include all np 
states up to n < Umax = 20. We find that at the location where Vf/^(r) = E{r) + 1 is a minimum in r, E{r) is also a 
minimum with respect to variations in 77 at ry = 0. We conclude that the local potential minimum of Vf/jj(r) at r ~ 6 
a.u. is stable against perturbations from e+-e~ correlations. 



VI. DISCUSSIONS AND CONCLUSIONS 



The four-body system of (e+e~pp) contains a large degrees of freedom as the particles can appear in many different 

combinations and correlations, exemplified by the large number of parameters in an elaborate variational calculation. 
To be able to classify the states of the greater (e+e^pp) system into simple families, if at all possible, will bring us to a 
better understanding of the complexity of the spectrum and the interactions that are associated with the complicated 
four-body problem. 

In terms of the correlation between the 6+ and the e~, one can perhaps roughly classify the greater {e^e^pp) system 
as belonging to the H-H family without e+-e^ correlations and the Pspp family with strong e+-e^ correlations. As 
one notes from energy considerations in Fig. 4, the lowest-energy state of the {e'^e~pp) system is dominated by a 
state of the H-H family at large interatomic separations but at short distances the system is dominated by a state of 
the Pspp family. The lowest-energy state therefore changes its character as the interatomic separation decreases. 

In the collision of the H and the H, the four-body system is prepared in the entrance channel of the H-H family as the 
composite particles approach each other. The system may follow the configuration of the lowest energy adiabatically 
in a slow collision, or it may follow the diabatic configuration of the H-H entrance channel when the collision is fast or 
when the intrinsic differences in the characteristics of the two configurations hindered their transition. In general, for 
the collision of arbitrary energies, it is necessary to couple the different channels explicitly and to follow the transition 
of the channels as the composite particles approach each other. We are therefore motivated to examine the interaction 
potential between the H and the H as composite particles, for which the e~ and the e+ remain uncorrelated, orbiting 
around the p and p respectively. 

We use the tool of the perturbation theory to examine this interaction potential, with the states of H and H as 
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unperturbed basis states. The perturbation theory can be a useful concept because the residue interaction Vj contains 
many cancellations and is small in magnitude compared to the value of typical energy denominators. The lowest-order 
perturbation theory gives an interaction potential Vdir{r) that is positive at large separations with a barrier located 
approximately between 1.4 to 5.0 a.u. and a precipitous drop at r < 1.4 a.u. Our calculations have been carried out 
up to the next-to-leading order, by including virtual excitations of bound and continuum / = 1 states of the H and 
H atoms. The next-to-leading order leads to the polarization potential that is always attractive. The bound-bound 
contribution to the polarization potential has been obtained by using the Fourier transform method involving the 
virtual excitation from the Is states to all np states with n < n,„ax = 20. The bound-continuum, and the continuum- 
continuum contributions have been obtained numerically by discretizing the continuum. The polarization potential 
results agree with previous calculations of the Cg coefficient at large separations. The variation of the interatomic 
separation then allows us to obtain the polarization potential Vpoi{r) at different separations. 

We find that the interaction potential between the composite H and the H possesses peculiar features. At inter- 
mediate distances, the total interaction potential Vfjf[{r), whic;h is the sum of the direct potential Vdir{r) and the 
polarization potential Vpoi{r), becomes negative. At shorter distances, the total potential shows a minimum at r ~ 6 
a.u., a barrier rising at r ~ 5 a.u., and a broad barrier peak at r ~ 1.9 a.u., as shown in Fig. 4. 

Our investigation on the stability against variations in the e+-e^ correlation in Section V indicates that the com- 
posite system is stable against variations of the e+-e~ correlation for r > 3 a.u. but is unstable for r < 3 a.u. Thus, 
the state with e+-e~ correlations is the state of lowest energy at distances r < 3 a.u., but the composite H-H state 
may remain stable without substantial e+-e~ correlations as an excited system for r > 3 a.u. 

The shallow local minimum located around interatomic separations of r ~ 6 a.u. has a depth of 6 x 10^^ a.u. One 
can use the interaction potential Vuni''') and solve the Schrodigner equation (12) for possible two-body H-H states. 
Even though the potential minimum is shallow, a simple WKB calculation indicates that the potential is deep enough 
to hold a bound state because of the large reduced mass of the hydrogen and the antihydrogen. The two-body H-H 
state, if it exists, will not be the ground state, as it is specialized to the configuration of the H-H family and lies 
above the state of the lowest energy obtained in a variation calculations. 

It is necessary in the future to include a greater number of the degrees of freedom, to take into account additional 
types of polarization of the two-body states such as the virtual excitation of the d, /,... states. As the residue interac- 
tion is dominated by the dipole-dipole interaction, the higher multipole interactions are expected to be diminishing in 
their importance and will not modify greatly the gross features of the interaction potential. Because the polarization 
potential is always negative, the addition of these multipole interactions will probably only deepen slightly the local 
potential minimum at r ~ 6 a.u. 

Returning to the question of channel coupling and reaction processes, we wish to emphasize that we have obtained 
only the interaction potential V^h for the channel of the H and the H as composite objects. We have in effect 
constructed the interaction potential only for one of the diabatic basis configurations. Interaction potentials for other 
basis configurations can be obtained in different analyses and a full pic;ture of adiabatic or diabatic motion for reactions 
at an arbitrary energy needs to be investigated by a couple-channel analysis involving different diabatic channels with 
different interaction potentials and transition matrix elements. 

In conclusion, we have examined the interaction potential of the H-H system as composite objects using the second- 
order perturbation theory with mutual polarizing excitations of the atoms. We found that the interaction potential 
has a peculiar local minimum at r ~ 6 a.u. and a barrier rising at r < 5 a.u. The potential energy minimum is found 
to be stable against perturbations of the e+-e^ correlation. Further definitive studies on the nature of the interaction 
potential energy around this energy minimum and the interplay between the H-H state and other states of different 
configurations will be of great experimental and theoretical interests. 
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APPENDIX A: FOURIER TRANSFORM OF THE TRANSITION DENSITY p^o^(r) 

We here derive the expression for the Fourier transform of the transition density, p^^{p), which involves the 
product of Is and np wave functions. The hydrogen wave function is given by 



i>mr.ir) = TV. ^fj^^) oxp J eV-r ) Y^r), (Al) 

where Nni is the normahzation constant, is the spherical harmonic and L^\x) is the Laguerre polynomial. 

Setting Z = 1, I =1 and ao — I, we have 

Pxo^i'^) = 4>ooo{r)(pnimir) 

= NooNni e-'-("+i)/"L3_2 j^^^ Foo(r)yi„(f). (A2) 

Prom Ref. [15], we have 

m = E f ^' + " ) - 0f Lf\{x). (A3) 



fe=0 



Now by letting r' = {n + l)r and writing 



we use Eq. (A3) with j = n — 2, a = 3, and /3 = l/(n + 1) , and we obtain 

n-2 



Wr') = E("J')/?'""'"'^(l-/?)' (v) ^'-^-'^ (?) • (A5) 



Therefore, the transition density Pxo^ {r') becomes 

Pxo'^ir') = ^^U,ir')Yoo{f)YUf). (A6) 
One can Fourier transform Pxq^ {f') by 

Pxo iP) = J—^^ 9x6 (P)' (A7) 

where 

ato'iP^ = j^Jdr'exp{ip.r'/{n+l)}U{r')Yim{f) 

= / rfr-'e'P'-^'Ki,/„ifc(r')yi^(f)/<i„ (A8) 

with p' = p/{n + 1). By recognizing the latter integrand is identical to the "p wave functions", without laborious 
mathematics and algebra, one can make use of the formula for the momentum space "radial" wave function of hydrogen 
atom given in Bransden and Joachain [18] and cast the Fourier transform of the above Px(,^{r') into the results of 
Eqs. (23)-(25). 



[1] M. Amorctti et al, Nature, 419, 456 (2002); Phys. Lett. B, 578, 23 (2004). 

[2] G. Gabrielse et al, Phys. Rev. Lett., 89, 213401, 222401 (2002); Adv. Atom Mol. Opt. Phys., 50, 155 (2005). 
[3] C. Y. Wong, A. K. Kerman, G. R. Satchler, and A. D. MacKellar, Phys. Rev. C29, 574 (1984). 



13 



[4] B. Zygclmau, A. Saenz, P. Froclich, and S. JonscU, Phys. Rev. A 69, 042715 (2004). 
[5] W. Kolos et al, Phys. Rev. A 11, 1792 (1975). 

[6] E. A. G. Armour, J. M. Carr, and V. Zcnian, J. Phys. B 31 L679 (1998); E. A. G. Armour, and V. Zeman, Int. J. Quantum 
Chem. 74 645 (1999). 

[7] P. Froelich, S. Jonsell, A. Saenz, B. Zylgelman, and A. Dalgarno, Phys. Rev. Lett. 84, 4577 (2000) 
[8] K. Strasburger, J. Phys. B 35, L435 (2002). 

[9] L. Labzowsky, V. Sharipov, A. Prozorov, G. Plunicn, and G. Soff, Phys. Rev. A 72, 022513, (2005). 
[10] V. Sharipov et al, Phys. Rev. A 73, 052503 (2006); Phys. Rev. Lett 97, 103005 (2006). 
[11] C. Y. Wong, Phys. Rev. C 69, 055202 (2004). 
[12] D. L. Morgan and V. W. Hughes, Phys. Rev. A7, 1811 (1973). 
[13] C. Y. Wong and H. W. Crater, Phys. Rev. C63, 044907 (2001). 

[14] G. R. Satchler and W. G. Love, Phys. Rep. 55, 183 (1979).; F. Petrovich, Nucl. Phys. A251, 143 (1975). 

[15] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover Pubhcations, Inc, New York) p785, 
eq. (22.12.7). 

[16] R. Eisenschitz and F. London, Zeitschrift fiir Physik. 50, 24 (1928); F. London, Nature, 17, 516 (1929). 
[17] K. Strasburger, J. Phys. B 37, 4483 (2004). 

[18] B. H. Bransden and G. J. Joachain, Physics of Atoms and Molecules (Longman Scientific and Technical, Longman Group 
UK Limited, 1992 ) p626, eq. (A5.34). 



